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NOMENCLATURE 


Symbol 


Description 



Gust spectrum multiplicative constant 
Vertical gust forcing coefficient vector 

Flap forcing coefficient vector 
Horizontal gust forcing coefficient vector 

Aspect ratio 

Span of one auxiliary wing (m) 

Mean aerodynamic chord of wing (m) 

Chord of flap (m) 

Lift coefficient at equilibrium 

Lift coefficient of horizontal stabilizer , based 
on wing area 

Lift curve slope of basic aircraft 
Lift curve slope of wing 

Lift curve slope of horizontal stabilizer, based 
on wing area 

Auxiliary wing lift curve slope 

Lift coefficent derivative due to downwash lag for 
basic aircraft 

Lift coefficient derivative with flap deflection 

Tail off lift coefficient derivative with flap 
deflection 

Lift coefficient derivative due to flap-induced 
downwash lag 

Lift coefficient derivative with normalized pitch- 
ing rate 


V 


Lift coefficient derivative with normalized horizontal 
gust acceleration for basic aircraft 

Lift coefficient derivative with elevator deflection 


Auxiliary wing lift coefficient derivative with auxili- 
ary wing tab deflection 

Pitching moment coefficient with angle of 
attack for basic aircraft 

Tail-off pitching moment coefficient derivative with 
angle of attack for basic aircraft, measured at center 
of gravity 

Pitching moment coefficient derivative due to downwash 
lag for basic aircraft 

Pitching moment coefficient derivative with flap de- 
flection, measured at center of gravity 

Tail-off pitching moment coefficient derivative with 
flap deflection, measured at center of gravity 

Tail-off pitching moment coefficient derivative with 
symmetric aileron deflection, measured at center of 
gravity 

Pitching moment coefficient derivative with elevator 
deflection 

Damping in pitch coefficient 

Pitching moment coefficient derivative with normalized 
horizontal gust acceleration for basic aircraft 

Pitching moment coefficient derivative due to flap 
downwash lag 

Longitudinal force coefficient derivative with angle of 
attack for basic aircraft 

Longitudinal force coefficient derivative with normal- 
ized speed perturbation 

Longitudinal force coefficient derivative with flap 
deflection 

Longitudinal force coefficient derivative with elevator 
deflection 

Flap hinge restraining moment coefficient required for 
horizontal gust alleviation 
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Flap hinge moment coefficient derivative with angle 
of attack. Includes effect of auxiliary wing. 

Hinge moment coefficient derivative with angle of attack 
for flap alone 

Flap hinge moment coefficient derivative with flap de- 
flection 

Flap hinge moment coefficient derivative with flap 
tab deflection. 

Hinge moment damping coefficient due to flap alone 

Flap hinge moment coefficient derivative with pitching 
jrate 

Drag coefficient derivative with flap deflection 

Thrust coefficient derivative with speed perturbation 

Measure of high frequency pitching rate response 
Total derivative with respect to x 
Aileron-flap coupling ratio for analytic model 
Aileron-flap coupling ratio for test data model 
Applied flap hinge moment for equilibrium(N-m) 

Flap hinge restraining moment required for horizon- 
tal gust alleviation (N-m) 

Flap hinge moment derivative with flap def lection (N-nv'r ad) 

Normalized moment of inertia about pitch axis 

Normalized moment of inertia of auxiliary wing about 
its longitudinal hinge line 

Reference normalized auxiliary wing moment of inertia 
Normalized flap moment of inertia 

Normalized flap-auxiliary wing system moment of inertia 

Incidence of auxiliary wing 

Vertical gust alleviation gain constant 

Vertical gust alleviation gain constant if no auxiliary 
wing is employed 
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Horizontal gust alleviation gain constant 
Arbitrary size of auxiliary wing 

Scale length of atmospheric turbulence spectrum (m) 

Distance between wing and stabilizer aerodynamic centers 
in mean aerodynamic chords 

Mass of aircraft (Kg) 

Zero degree coefficient of characteristic equation 

2 

Normalized pitching rate; also dynamic pressure (N/m ) 
Radius of spring linkage from hinge line (m) 

Wing area (m^ ) 

2 

Horizontal stabilizer area (m ) 

2 

Auxiliary wing area (one wing) (m ) 

2 

Flap area (one flap) (m ) 

Static margin 

Thickness of auxiliary wing components (m) 

Normalized speed perturbation 
Normalized horizontal gust velocity 

Aircraft velocity (m/sec) 

Aircraft weight (N) 

Angle of attack perturbation 

Angle of attack perturbation due to vertical gust 

Angle between flap linkage radius and flap datum 
Auxiliary wing to flap gearing ratio, 36 /36 

3.V7 

Flight path angle 

Normalized rate of change of flight path angle 
Flap deflection 

Auxiliary wing deflection about longitudinal hinge line 
Elevator deflection 

Flap or auxiliary wing tab deflection 
Load factor (g) 



3e/3a Downwash derivative with angle of attack for basic 

aircraft 

3e/36 Downwash derivative with flap deflection 

C Damping ratio 

0 Pitch attitude 

VI Normalized mass, 2m/pSc 

p Air density (Kg/m^) 

a Root mean square gust velocity (m/sec) 

$ Power spectral density of pitching rate 

$ Power spectral density of vertical gusts 

Wg 

Power spectral density of load factor 

Spacial frequency (rad/m) 

(D Temporal frequency (rad/sec) 

Normalized natural frequency 

0)^ Normalized natural frequency of free-floating flap 

Normalized phugoid natural frequency 
P 

w Normalized short-period natural frequency 

sp 


Superscripts 

' Pertains to alleviated aircraft 

Subscripts 

f Applies to flap 

aw Applies to aiixiliary wing 
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I . INTRODUCTION 


The study has been undertaken to determine the feasibility 
of improving ride comfort in a light aircraft in cruise without 
having to resort to the use of relatively expensive servomecha- 
nisms. The system hopefully will be attractive to industry to 
offer as a low cost option and, perhaps, even as a retrofit pack- 
age to present aircraft. The major design considerations for 
different aircraft center on modifying the flap design and as- 
certaining flap hinge moment characteristics, from which the size 
of the auxiliary wings can be prescribed. 

When attempting gust alleviation, there are several basic 
problems encountered which lead to distinct design and perform- 
ance compromizes. Alleviation implies maintaining a constant 
load factor which is normally influenced by vertical and longi- 
tudinal gust components. At low cruise lift coefficients, loads 
due to vertical gusts predominate by a significant margin and 
horizontal gusts can be neglected. This would not be the case, 
for instance, of a STOL aircraft taking off or landing. Com- 
pensating both gust components introduces some interesting sta- 
bility problems which are discussed briefly later. Alleviating 
vertical gusts corresponds to maintaining a constant lift coef- 
ficient. Consequently, the lift curve slope is zero and the 
damping in pitch is reduced. 

The most practical way to manipulate the lift coefficient 
is with trailing edge flaps. Although their linear range of 
operation is limited, they nevertheless can compensate gusts up 
to about 3 m/sec which covers most of the turbulence commonly 
encountered. However, unfavorable pitching moment and downwash 
characteristics are introduced by these devices. The increased 
camber produces a nose-down pitching moment and the modified 
downwash pattern has a strong influence on the downwash deriva- 
tive with respect to angle of attack at the tail. Although for 
normal inboard flaps the latter effect somewhat offsets the pitch- 


ing moment due to camber, the net result is a sharp reduction in 
pitch stiffness ) . Downwash lag also is affected, and results 

in lower damping in pitch and amplified pitching response to gusts. 
One cannot design a flap only system which will maintain a constant 
load factor and zero pitching rate response while providing 
adequate stability. 

There is one more difficulty associated with the dynamic 
behavior of gust-alleviated aircraft which, fortunately, can be 
circumvented by proper design. Because elevator control is effec- 
tive through the influence of the pitching moment on angle of 
attack, and because the lift curve slope of the alleviated air- 
craft is zero, the aircraft is not controllable in the normal 
sense. The flight path angle will not respond satisfactorily 
to elevator movement. Some means of overriding the system must 
be provided to restore maneuverability- 

The basic concepts of gust alleviation are not new. In 
1938 Rene Hirsch^ of France conducted model tests of a flap-type 
alleviation system which employed the horizontal stabilizer as 
an angle of attack sensor. These surfaces were articulated on a 
longitudinal axis and coupled to the flaps. Hirsch was the first 
to recognize the need to concentrate on long wavelength absorb- 
tion and was not concerned with the sensor lag of this configura- 
tion. His interest continued after World War II when he built a 
small twin-engined aircraft using these concepts. He attempted 
to absorb both vertical and horizontal gusts, but found it neces- 
sary to avoid the zero root attendant this system by employing 
altitude rate coupling through a variometer to the elevator. 

His choice of mechanization was rather ingeneous in that it auto- 
matically bypassed the alleviation system and provided the maneu- 
verability associated with a conventional aircraft. 

In the early 1950' s, Phillips and Kraft^ investigated the 
effects on the short period dynamics of introducing a vane-type 
angle of attack sensor mounted on the nose of a twin-engined air- 
craft and driving the flaps and elevator through a hydraulic 
servo. They developed the relationships identifying the flap 
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characteristics required to provide vertical gust alleviation 
and zero pitching response and recognized the degeneration of 
the short period roots to a real pair, one being zero and the 
other having a very short time constant. In addition they 
pointed out the difficulties of designing a flap system to 
achieve the specified conditions and suggested the possibility 
of using a segmented flap to help meet the downwash requirements 
The concept was applied to a full scale aircraft^ which demon- 
strated the feasibility of adjusting the flap downwash to im- 
prove the pitching response of the alleviated aircraft. Maneu- 
verability not only was recovered, but was improved by coupling 
the stick to the flap actuator, thereby overriding the allevia- 
tion system and providing direct lift control. In a recent re- 
port by Phillips^, the early analysis was extended to include 
horizontal gust alleviation and its effect on the phugoid mode. 
He confirmed the stability problems introduced by such a system 
as first recognized by Hirsch. He further examined the horizon- 
tal accelerations resulting from variations in drag due to the 
flap as it responds to the gusts. In a numerical example of a 
STOL aircraft, these accelerations amounted to as much as sixty 
percent of the values for the vertical accelerations which would 
be experienced by the unalleviated aircraft. A method of coup- 
ling spoilers to the flaps was suggested as a means for re- 
ducing the net drag due to flap deflection, thereby diminishing 
the source of these accelerations. 

The present system under consideration is forced to com- 
promize the optimum alleviation conditions not only to maintain 
stability and handling parameters, but also to work within the 
framework of system objectives, namely, to make use of existing 
control surface regions. The flaps must be redesigned to permit 
upward deflection; a plain, mass-balanced flap with hingeline 
at twenty percent chord is substituted. Furthermore, while anal 
ysis for both vertical and horizontal gust alleviation is pre- 
sented, only the former is proposed for a conventional light air 
craft. 
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The system is analyzed for the Cessna 172 aircraft. The 
angle of attack sensors are chosen to be two auxiliary wings, 
large enough to drive the flaps directly and located ahead of 
the cargo compartment door so as not to impede entry to the cabin. 
This concept was originated by W. H. Phillips as an alternative 
to the Hirsch method of using the empennage as the gust sensor. 
These wings are mass balanced and are articulated in roll to 
deflect the flaps and in pitch for coupling to the elevator in 
order to regain basic aircraft maneuverability. Elevator coupling 
with the flaps is not considered necessary to reduce pitching 
motions (although it could be required for other aircraft) ; com- 
plication of the control stick - elevator interface is thereby 
avoided. 



Figure 1. Auxiliary Wing Installation 
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The bulk of the analysis presented here is predicfited on sev- 
eral fundamental assumptions regarding the behavior of the air- 
craft. These are: 

1. The aircraft is a rigid body. 

2. Unsteady aerodynamic effects can be neglected. 

3. The flap operates only in its linear lift curve region. 

4. The equations of motion are linear in the perturbation 
variables. 

5. The auxiliary wing-flap system responds instantaneously 
compared to the rigid body motions. 

The first two are justified by considering the spectra of gusts; 
the power is negligible at frequencies much higher than that of 
the short-period mode. The third and fourth merely impose limits 
on the region of validity for the analysis. The fifth is made 
to simplify the analysis and to identify the capability and 
limitations of an idealized flap-auxiliary wing system whose 
dynamic characteristics do not influence the gust response of the 
aircraft. Actual responses for a wide range of flap system dynamic 
parameters are presented to call attention to the importance of 
carefully selecting system components to achieve satisfactory 
response characteristics. 

If the flap characteristics are considered to be linear 
over a range of fifteen degrees, they can nullify an incremental 
lift corresponding to about three degrees of change in angle of 
attack. At cruise speed, this value translates into a vertical 
gust of about three meters per second. 
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II. SYSTEM DESIGN 


The auxiliary wings must be large enough to provide the 
power necessary to drive the flaps but care must be exercised 
not to overalleviate the gusts since this would lead to insta- 
bility. While it is recognized that theoretical estimates of 
flap hinge moment parameters involve substantial uncertainty and 
whenever possible , should be replaced by experimental data from 
full-scale wind tunnel tests or actual flight information, they 
nevertheless permit sizing estimates. 

The steady-state response of the flap to gusts can be des- 
cribed by a static hinge moment equation in which moments due to 
gusts are balanced by that developed by a flap deflection. 


The last term results from spring loading the flaps which creates 
a torque independent of flap deflection and dynamic pressure. 

The required hinge moment coefficient derivative with angle of 
attack can be formulated from the first two terms of equation (1) 
and equation (ElO) which defines the relation between flap de- 
flection and gust induced angle of attack for optimal gust alle- 
viation. 


'H 


= C 


a 




(2) 


If a pushrod is used it is attached to the auxiliary wing quarter 
chord line at the point midway along the span to neutralize the 
aerodynamic loads sustained by the hinge supports. For small 
amplitude motions the auxiliary wing lift is independent of flap- 
ping angle. Consequently the only flap parameter to be substan- 
tially modified by the coupling is the hinge moment sensitivity 
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to angle of attack, 
meters 



In terms of flap and auxiliary wing para- 


= C 


H 


a 


, 36 eb) 

1 aw aw _aw 

2 36 (Sc)j 


(3) 


In view of its small size the auxiliary wing contribution to C_ 

Ij 

is assumed negligible. ^ 

Equations (1) and (2) can be combined and solved for the 
span of the auxiliary wing, if the aspect ratio and airfoil sec- 
tion characteristics are known. An aspect ratio of 2.0 and an 
NACA 0009 airfoil yield a lift curve slope of 2.54. (Ref. 9) 

The gearing ratio between the flap and auxiliary wing was select- 
ed to be unity. The resulting wing span is 1.13 m. A larger 
gearing ratio would lead to a smaller auxiliary wing size but 
larger deflections (6 ) . 

Horizontal gust alleviation is provided by loading the 
flaps to achieve the sensitivity specified in Equation (E-20) . 





(4) 


If a spring is used to realize this moment, it effectively must 
have zero stiffness to avoid changes in torque due to flap deflec- 
tion, yet must be able to be set to the torque level specified 
by the weight of the aircraft, which varies during flight. 

In equilibriiam flight, the flap deflection can be trimmed 
to zero by adjusting the incidence of the auxiliary wing to counter- 
act the normal floating tendency of the flap. Thus, engaging the 
gust alleviation system would not affect the aircraft trim condi- 
tion. If horizontal gust alleviation is not employed, this inci- 
dence setting is independent of airspeed. 
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III. MODIFIED AIRCRAFT DYNAMICS 


A relatively simple set of equations of motion can be used 
to investigate the dynamic behavior of a gust-alleviated aircraft. 
The usual longitudinal equations of motion are augmented with the 
flap-auxiliary wing system hinge moment equation. 



To simplify the analysis and to focus attention on the modifica- 
tion of the rigid-body modes introduced by the gusts alleviation 
system, the flaps are assumed to operate fast compared with the 
rigid-body motions and the flap equation can be reduced to a 
static balance similar to that given in equation (1) but in- 
cluding the perturbation variables a and u . With this simpli- 
fied equation, the conditions for gust alleviation (equations E-10 
and E-20) are introduced. It is convenient to introduce gains 
representing degrees of alleviation, which are useful in examining 
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the effects of the system on the aircraft dynamics, 
equation becomes ^ 

S ' -KaC^<“+“g'-Kv^<»+“g) 


The flap 
( 6 ) 


The common assumptions associated with small perturbation 

techniques for conventional aircraft apply. ^ In addition, for 

the low Mach number regime of the aircraft, C is assumed zero. 

• ^u 

Lift and pitching moment derivatives with u^ are included as 
forcing terms to examine the response of the aircraft at fre- 
quencies just above the short period mode. 


Ci -2ys 
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a a q 

P 

1 

u 

2C ^ 

-K + C s 

vC^ m,. m* 

L- 0 u 
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(7) 


The primed quantities are interpreted as follows: 

■ 

C = C -K 

x^ aC- X. 
a a Ljc o 


( 8 ) 


C = C -K 
x„ x„ VC- Xc 

UL U Xj o O 


(9) 
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For the lift and pitching moment terms, L and m are s\abstitut- 

ed for X in these expressions. The rate-dependent derivatives 

are similar except that the last parameter in each expression 

• * . 
is a function of 6 . 

To comprehend the consequences of gust alleviation, these 
equations should be examined without imposing all of the optimal 
conditions simultaneously. The lift curve slope is reduced to 
zero for K =1, thereby eliminating the constant term in the 
angle of attack perturbation of lift. The corresponding velocity 
perturbation term also is eliminated if K^=l. The combination 
of the two conditions leads to a zero root for the characteristic 
equation when the nominal flight path angle is zero. 

The static stability of the aircraft is determined by the 
sign of the constant term P(0) in the expansion of the character- 
istic determinant. 


P(0) = 



DC 


M 


Da 




(1-K^)] 


(10) 


It is thus apparent that the fully alleviated aircraft is mar- 
ginally stable since both of the alleviation constants are unity 
and P(0)=0. If full horizontal gust alleviation is only partially 
implemented (K_=l, K <1) 

Si V 

9 DC 

P(0)= -2C^ (1-K^) (11) 

SO that in this case, static stability requires negative total 
pitch stiffness. 


DC 


M 


Da 


<0 


( 12 ) 
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The following table illustrates the possible range of this ■ 
parameter for the Cessna 172. 


Configuration DC_/Da 


flap -0.15 

aileron .93 

flap-aileron .40 


Thus only the flap configuration can provide static stability/ 
for the stability derivatives assumed. Test data recently supplied 
by NASA indicate the flap pitching moment derivative is smaller 
than estimated by this analysis, thus expanding the possibilities 
of stable operation. 


If the flap pitching moment characteristics are developed. 


C' can be expressed by 


m 


= ^m -^ac, 
a a L 


a 


[C 


m. 


-C 


m 


w 


a. 


3c . 

3^J 


( 13 ) 


For a positive flap downwash derivative and center of gravity 

behind the wing aerodynamic center, the large negative pitching 

moment due to increased camber is somewhat offset, but never is 

eliminated completely. If the flaps were designed to produce a 

negative downwash derivative, as specified in Appendix E, 

the aircraft would be unstable unless the elevator was coupled 

to the flap to provide a compensating pitching moment. The 

independence of C! on the center of gravity location is developed 

easily. The trimmed, alleviated aircraft exhibits no pitching 

moment about the center of gravity. If the center of gravity, is 

moved, the resulting moment is equal to the lift multiplied by 

the distance moved. The corresponding change in C' is this 

““a 

distance multiplied by C' . However, the latter is zero, 
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resulting in the independence cited. Hence, the static stability 
of the vertical gust alleviated aircraft cannot be improved 
by moving the center of gravity forward as is possible with a 
conventional configuration. A similar independence on center 
of gravity location is found for when horizontal gust 

alleviation is mechanized. 

C^, depends on the downwash lag and is influenced strongly 

by 6e/3^*. For inboard flaps the flap downwash derivative is large 

and positive, and reduces the short period damping. For the 

optimal case, it is negative and is increased to the value 

of C . 
ni 


a 


As more of the optimal alleviation conditions are imposed, 
more of the coefficients of the gust forcing terms are eliminated. 
When they all are satisfied, no forcing terms remain. As this 
process proceeds, assessment of the dynamic behavior of the air- 
craft becomes complex and at this point a numerical example is 
introduced to show the results of a practical application. The 
locus of roots for a Cessna 172 is shown in Figure 2 where K 

cl 

is varied from 0 to 1.5 and is zero or one. The short period 
mode is independent of and decreases in frequency and damp- 
ing with increasing K^. If loci for a variety of static margins 
were superimposed, they all would intersect at the point where 
K «1, since C' =0 and C is’ independent of static margin. The 
stability of the alleviated aircraft is demonstrated by the 
phugoid mode loci. For K^=0, the aircraft is tolerant of minor 
over-alleviation and yields roots which can be accomodated 
by the pilot. However, for K^=l, the system is extremely sensi- 
tive and diverges rapidly if over-alleviated by even the slight- 
est amount. A companion locus, not presented, shows somewhat 
less sensitivity to over-alleviation of horizontal gusts when 
K =1. The loci for the aileron-configured aircraft (Figure 3) 
where outboard flaps reverse the sign of the downwash derivative. 
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Stotic margin «0 


Flop oonfigurotien 



Stotic margin ■ 0 Flop configorotion 



Figure 2 Root Locus of Flap-Configured Aircraft as a 
Function of Vertical Gust Alleviation Gain 
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Rod/MC 


Stotic margin ».05 


Flop configuration 



Stofic margin ■ .05 Rop conf iguration 



Figure 2 continued 


Rod/MC 









Static morgin >0 


Ailofon configurotlon 



Stotic margin * O Aileron configuration 



Figure 3 Root Locus of Aileron-Configured Aircraft as a 
Function of Vertical Gust Alleviation Gain 
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static margin* .1 


Aiitron configuration 



Static margin ■ .1 Ailoren configuration 



Rad/stc 

Figure 3 continued 
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show that very little vertical gust alleviation 
can be mechanized before the onset of instability. A similar, 
less stringent condition is imposed on an aircraft fitted with 
full-span flaps. (Figure 4) . Recent test data generated by 
NASA indicates the tail-off flap pitching moment derivative with 
deflection has about half the computed value. An analysis of 
the short period mode is presented in Chapter VII and shows 
both the aileron and full-span flap configurations to be 
stable for K =1, K =0. This analysis clearly emphasizes the 
need for accurate parameter estimation before undertaking 
alleviation system design. 

Because the stability of the aircraft is very sensitive to 
over-alleviation of either vertical or horizontal gusts if the 
other is fully implemented, it is appropriate to carefully assess 
the relative influence of these two components on ride comfort. 

For equal gust amplitudes the ratio of load factors due to ver- 

ical and horizontal gusts is C, /2C^ which suggests very large 

ct 

lift coefficients must be realized before alleviation of horizon- 
tal gusts becomes important. It is for this reason that only 
vertical components are considered for light aircraft in cruise. 

For the Cessna 172, using inboard flaps as gust alleviators, 
C' remains negative. As shown in the next section, the pitching 
response to gusts is increased. This problem could be moderated 
by proportionately coupling the ailerons to the flaps to the 
extent that the aircraft remains stable. This technique provides 
the designer with a means of controlling the flap downwash at 
the tail. However, it introduces the additional complexity of 
developing linkages to permit both symmetric and antisymmetric 
deflections. Another degree of control can be exercised by 
coupling the elevator to the flaps as proposed by Phillips. In 
the case of all mechanical linkages, though, some thought must 
be given to the cockpit interface to prevent undue loading of 
the stick. 
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Rod/$»c Rod/MC 




Rod /sec 

Figure 4 continued 
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IV. RESPONSES TO GUSTS 


The Cessna 172 application is used to examine the frequency 
responses of basic and alleviated aircraft. It is important to 
recognize that these are not responses of an optimally alleviated 
aircraft; rather they are the result of working within the confi- 
guration constraints mentioned previously. However, the response 
levels and ramp characteristics cited parametrically on the 
graphs are general and can represent the optimal case. These re- 
lations can be developed by making use of the short period and 
phugoid approximations^ and assuming the pitching rate, angle of 
attack rate, and longitudinal gust rate derivatives result only 
from tail loads. The first four figures represent an idealized 
configuration, with instantaneously responding flap; the last seg- 
ment of this chapter presents curves for a wide range of flap 
system characteristics. 

The load factor response (Figure 5) of the basic aircraft 
is greatest at frequencies above the short period. There are 
two real zeros, one near this frequency and the other above 100 
rad/sec which produce a flat response in this band. The inertia 
of the aircraft dominates and prevents any significant buildup 
of vertical speed to offset the gust-induced angle of attack and 
complete loading, proportional to dL/da, is realized. At lower 
frequencies, aircraft motion does offset the gusts, except at 
the phugoid frequency where the very low damping produces a sharp 
peak in the response. At the lowest frequencies, the response 
is determined by the convection condition, whereby the aircraft 
moves with the air mass and the total angle of attack perturba- 
tion is zero, a+a =0. 

For the alleviated aircraft, the zeros transform to an un- 
damped pair just above the short period frequency, which intro- 
duces an upward slope in the response at high frequencies. This 
represents tail loading only and results from not satisfying the 
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Figure 5 Frequency Response of Load Factor to Vertical Gusts 



Figure 6 Frequency Response of Pitching Motion to Vertical Gusts 
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downwash criterion. The short period is longer because the lift 
curve slope is zero and C is diminished. At low frequencies, 

mot 

the phugoid persists when only vertical alleviation is implement- 
ed, but disappears with full alleviation. For the Cessna 172, 
gust alleviation is quite good in the region of interest, namely 
between the phugoid frequency and ten rad/sec. Above that the 
analysis is of limited value because unsteady lift effects, flap 
dynamics, and bending modes have been neglected. These frequen- 
cies are unimportant with regard to gust loads because the gust 
spectrum is very weak in this region. The effects on the air- 
craft are further diminished at the highest frequencies by the 
isotropic nature of these gusts which affords some spanwise 
averaging of the induced loads. 

The significant pitching behavior (Figure 6) can be de- 
duced from the response at very high frequencies. Inertia pre- 
vents any angular perturbation of the aircraft and the pitching 
moment results only from gust rate dependent tail loads. Con- 
sequently, the pitching rate is in phase with the gust angle of 
attack. For the alleviated Cessna 172, the forcing term in- 
creases in magnitude as a result of the change in sign of C^J^^ , 
which is due to the flap downwash lag. It is difficult to “ 
estimate the effect of this pitching motion on ride comfort. 
Although it is three times greater than for the basic aircraft, 
gusts are weak at these high frequencies. A flight test evalu- 
ation can best determine its significance. 

The responses (Figures 7 and 8) to horizontal gusts for the 
basic aircraft are much smaller than their vertical gust counter 
parts except near the phugoid frequency. As would be expected, 
there are only minor differences between the responses for the 
vertical-gust-alleviated aircraft and those for the basic air- 
craft over the entire frequency range. These are due mainly to 
differences in mode shapes for the two aircraft. Implementation 
of horizontal gust alleviation eliminates the load factor forcing 
term dependent on u , leaving only the rate-dependent term which 
leads to the ramp response. The pitching rate response for this 
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Figure 7 Frequency Response of Load Factor to Horizontal Gusts 



Figure 8 Frequency Response of Pitching Motion to Horizontal Gusts 


34 






case is much greater" at high frequencies than for the basic air- 
craft since C , is nearly an order of magnitude larger than C 

"*6 % 
Despite this increase, the response is less than one third g 

the pitching rate response of the alleviated aircraft to vertical 
gusts. 

In order to redirect attention to the effect of flap dynamics 
and to investigate any possible interference with the short-period 
mode of the aircraft, the complete equations of motion were used 
to obtain the load factor and pitching responses to vertical 
gusts corresponding to a variety of flap-auxiliary wing system 
inertia and aerodynamic damping. The initial design, with a 
natural frequency of 13.2 rad/sec and damping ratio of .45 has 
a rather poor load factor attenuation (slightly in excess of 50%) 
at the basic aircraft short-period frequency. As the natural 
frequency of the system is increased while maintaining a constant 
damping coefficient (C„.=-2.5), for example by reducing the 
structural weight of the system without changing the auxiliary 
wing shape or size, very little improvement in load factor 
attenuation is achieved, (Fig 9). This is due to the fact that 
the damping ratio of the flap mode is inversely proportional to 
the natural frequency of the system; thus, as the inertia is 
reduced, critical damping is reached and the degenerate mode 
becomes essentially damping limited. Therefore, in order to 
shorten the flap-auxiliary wing response time, the natural 
frequency of the system should be increased and the aerodynamic 
damping reduced correspondingly. Examination of the relations 
developed in Appendix A shows that the inertia and damping 
contributed by geometrically similar auxiliary wings are inver- 
sely proportional to their size and size squared respectively, 
if the gearing ratio is chosen to achieve full vertical gust 
alleviation (K_=l, Fl =1). Thus in general terms, the larger 
the auxiliary wings, the better the load factor attenuation. 

The effect of the reduced damping is clearly seen on the 
load factor frequency response curves corresponding to a constant 
damping ratio ( 5 =. 45) equal to that calculated for the initial 
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flap-auxiliary wing design. Note that the zero inertia curve on 
(o)^=oo ) represents the idealized case of instantaneous flap deflec- 
tion. Even better results are obtained by going towards smaller 
damping ratios as seen on Fig 11 . Improved load factor attenuation 
in the short-period frequency range is achieved at the expense 
of an increased resonance at the flap system natural frequency 
where the gust energy is much smaller. A damping ratio of .2 
and flap-auxiliary wing natural frequency of 20 rad/sec seems to 
be a satisfactory compromize. Good overall load factor attenuation 
with nearly one decade at the basic aircraft short-period shows 
the real utility of the passive gust alleviation system. 

Furthermore these dynamic characteristics are believed to be 
achievable in practice by proper optimization of auxiliary 
wing weight size and shape. 

Since alleviation of normal accelerations is not optimal 
with the real flap system, the high frequency pitching response 
is not as severely aggravated as suggested by the simplified 
analysis. Fig 14 shows this response is only 1.6 times as large 
as that for the basic aircraft at the short-period for a flap 
system with m^=20 rad/sec. This should cause no hardship on 
passenger comfort. 
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Figure 10 Effect of Flap Dynamics on Load Factor Response 
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V. FLIGHT IN A TURBULENT ATMOSPHERE 


The one dimensional Von Karman spectrum is chosen to 
model the vertical component of high altitude turbulence. 
The vertical gust power spectral density is 


(S5) 


27T 


l+|(ajfn) ^ 
[i+(a;fn)^] 


( 14 ) 


a=l. 339 


a is the root mean square gust velocity and can be upwards of 
3 m/s in rough weather. The scale of the turbulence is large, 
typically a few thousand fleet and has a critical efflect on the 
relative levels ofl excitation of the two rigid body modes. As 
the scale of the turbulence is increased, the spectrum shifts 
towards the lower frequencies while increasing in level below 
the break frequency. The power spectral density of the normal 
acceleration is expressed in terms of the gust spectrum, 

g g 

The gust wave number n is related to the frequency 03 through: 

n = ^ (16) 

Similarly for the pitching velocity power spectrum, 

“I CO? (i(i))|^ $ (to) 

^ g g 

A convenient tool for easily evaluating the mean-square 
gust responses and for identifying the most significant fre- 
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quency ibands in these responses is a graph of (on a linear 

scale) versus frequency (on a logarithmic scale) . The resulting 
curve can be integrated linearly to yield the mean square value.' 

oo 

= jf$jj(a))doj = 2*Area‘lnlO (17) 

— 00 


By neglecting the contribution of gusts over 10 rad/sec an 
approximate value can be obtained for the basic aircraft's re- 
sponse to atmospheric turbulence with, for example, scale length 
of 305 meters. 


rms An 
a 


0.082 g-sec/m 


rms ( 
a 



0.65°/m 


(18) 


For the fully alleviated aircraft with instantaneously acting 
alleviation system (flap configuration, 15% 31), 


rms An 
a 


= 0 



2 . 3°/m 


(19) 


Power spectral density plots of load factor and pitching motion 
versus frequency of the type previously discussed (J2$ vs. Log w) 
are presented for several degrees of alleviation. The aircraft 
is flap- configured and responses for three values of basic air- 
plane static margins are examined. The scale of turbulence was 
chosen to be 305 m which sets the dominant gust wavelength in 
the phugoid range. 

As the gain of the vertical gust alleviation system is in- 
creased, the contribution of the high frequency gusts to the r.m.s. 
load factor is progressively reduced to zero while at the same 
time the pitching motion in the short period frequency band is 
considerably increased. In order to get any significant acceler- 
ation alleviation, K must be near unity. On the other hand the 
introduction of even a little vertical gust alleviation results 
in a sizeable increase in pitching response. 
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The horizontal gust alleviation influences only the low 
frequency response and essentially eliminates the resonant peak 
at the phugoid frequency. 
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Figure 15 Gust Spectral Responses 
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Figure 15 continued 
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VI. MANEUVERABILITY 


Analysis of the effects of the gust alleviation system 
on the maneuverability of the aircraft shows that the elevator 
control effectiveness is almost completely lost. From a prac- 
tical standpoint it is convenient to differentiate between short- 
term and long-term response to control deflections. The short- 
period mode is essentially a constant speed oscillation in which 
normal and pitching motions occur. Consequently, spring loading 
the flaps to attenuate normal accelerations produced by changes 
in dynamic pressure has no effect on the response to high fre- 
quency elevator deflections. Thus both the vertically -alleviated 
and the fully -alleviated aircraft have the same short term ma- 
neuverability. 

The zero value of the effective lift curve slope is respon- 
sible for the loss of short-term control of flight path angle. 
Because of the constant wing lift, the flight path angle is 
governed by the tail loading. When the stick is pulled back, the 
reduced lift on the tail results in the aircraft sinking slowly. 
The reduced pitch stiffness produces a considerable increase in 
the attitude and angle-of -attack sensitivities to elevator deflec- 
tions in the phugoid to short-period frequency band. For the 
Cessna 172 the angle of attack response is increased by a factor 
of eight so that even minor stick deflections might quickly lead 
to wing stall. 

The long term response of the aircraft is determined by the 
phugoid mode where horizontal gust alleviation has an important 
effect. The constant static lift condition, around which the fully 
alleviated aircraft is designed, causes the phugoid to degenerate 
into a pair of real roots, one of which is zero, so that the flight 
path diverges with a slightly untrimmed elevator. On the other 
hand, if just vertical gust alleviation is implemented the phugoid 
is affected only moderately. When the elevator angle is varying 
slowly, the angle-of -attack perturbation is no longer small and 
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is determined by the condition that pitch equilibrium is maintain- 
ed. The steady-state angle of attack is multiplied by a factor 
of six and the response is flat up to the reduced short-period 
natural frequency. Another negative aspect of the gust alle- . 
viation system is the loss of the airspeed trimming capability 
of the elevator which is due to the zero lift curve slope. The 
steady-state velocity perturbation changes sign and is reduced 
by a factor of twenty. 

The basic aircraft long term attitude response to a step 
of elevator occurs primarily at the phugoid frequency. A signi- 
ficant change in this response occurs when angle-of-attack gust 
alleviation is introduced. A pair of zeros in the attitude 
transfer function virtually eliminates the resonant peak of the 
phugoid so there is very little residual long-period oscillation 
in the transient response. The small pitch stiffness is respon- 
sible for the increased steady-state pitch angle. 

The maneuverability of the vertical-gust-alleviated aircraft 
is recovered without significant change in stability character- 
istics or sacrifice of riding qualities if the elevator and flaps 
are coupled to offset the hinge moments produced by changes in 
angle-of-attack so that the flaps remain in static equilibrium 
when the elevator is operated. This might be accomplished, for 
instance, by installing a tab on the auxiliary wing or on the 
flap itself, or by articulating the whole auxiliary wing about 
its quarter-chord line. The feasibility of each depends on the 
coupling and hinge moment requirements. With the addition of a 
tab, the hinge moment equation for the flap system becomes: 


'H, 


6 + C„ a + C 


'H 


a 


H, 


= 0 


( 20 ) 


The tab-elevator gearing 96^/36^ is determined by the con- 
dition that the flaps remain undeflected 




_gt_ 

6 


t 


0 


(21) 
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where (c*/5g)Q denotes the angle of attack transfer function of the 
basic aircraft. Since a mechanical linkage with a fixed gearing 
ratio is to be used, the control recovery condition cannot be sa- 
tisfied at all elevator frequencies. This is of no great impor- 
tance however, since (ct/6g)g essentially is flat up to the short 
period frequency. At low frequencies. 




1.52 


( 22 ) 


Between the phugoid and short period frequencies it is : 


in p 

, a ^ °e 2ySl 

'6 ^0 “~C 2y»I-C 
® "*0 "‘q 


- 1.11 


(23) 


Therefore if the gearing is chosen to recover low frequency con- 
trollability it also should be satisfactory at high frequencies. 
This is confirmed by the frequency responses (Figures 16 to 19) 
and by analog simulation of the aircraft response to step eleva- 
tor inputs (Figure 20) . Control of airspeed and flight path 
angle are virtually fully recovered. Angle of attack and pitch- 
ing response are increased at frequencies above the phugoid 
although this is somewhat offset by the reduced natural frequency 
of the short period motion. The main difference is to be found 
in the larger phugoid oscillation in angle of attack; however 
this is not perceived easily by the pilot. 

In order to keep tab deflection small while allowing for 
full elevator travel the gearing ratio should be kept reasonably 
close to unity. When Equation (2) is applied to the tab-eleva- 
tor gearing equation, the result is, 
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Figure 20 Time Responses to a Step Elevator Input 
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It is apparent that fitting tabs to the trailing edge of the 
flaps would require an excessively large gearing ratio between 
elevator and tabs for the present flap configuration which is 
not aerodynamically balanced. Inspection of the flap-auxiliary 
wing coupling relationships shows that fitting tabs to the 
auxiliary wings leads to a smaller gearing ratio, 

,aw 


c 

96^ ^L. H 


96 


t a , a . 

_aw M ^ 


a 


0 C -C^ 
6 ^ a a 


(25) 


When the flap contribution to the total angle of attack derivative, 

, is small compared to that of the auxiliary, wing, this relation- 
"a 

ship can be interpreted physically as meaning that stick deflec- 
tions have no effect on the auxiliary wing position so that the 
flaps remain undeflected. The expression can be rewritten. 



1.52 


(26) 


It is apparent from this relation that fitting a tab to the 
auxiliary wing might require too great a gearing ratio between 
tab and elevator to permit each to operate over the full extent 
of its linear range. The best solution is to tilt the entire 
auxiliary wing about its quarter chord line, thus changing its 
incidence, i^^f and reducing the corresponding gearing ratio to 
the minimum value of 


9i 


m, 


aw 


96. 


= 1.52 


m 


a 


(27) 


This effectively means tilting the auxiliary wing by an amount 
equal to the steady-state change in angle of attack due to ele- 
vator deflection. If is not neglected in this expression, 

“Q( 

the result is a slight tilting past the free stream angle to 
offset this torque. 
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VII. SENSITIVITY ANALYSIS 


Because of the large uncertainties in the estimation of the 
flap hinge moment parameters and the poor load factor attenuation 
of the under alleviated system, but more importantly, because of 
the generally unacceptable and possibly dangerous stability char- 
acteristics resulting from overalleviation, it is necessary to 
determine the relative sensitivities of the general alleviation 
constant K_ to the various system parameters and to use these to 
devise methods of safely implementing full alleviation. In addi- 
tion, knowledge of the behavior of the short period roots to un- 
certainties in flap pitching moment and downwash derivatives per- 
mits the designer to make full use of various control surface 
geometries and to determine the need to couple elevator with the 
flap to modify the pitching tendency. Recall that 



The augmented flap hinge moment derivative in the presence of the 
auxiliary wing is 


r = - iP (S b) aw aw 

2 I (S=), 


(29) 


where 6 


The relative sensitivity of K to C„ 

3. H t 


can 


be computed by taking the ratio of logarithmic partial derivatives 
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( 31 ) 


Similarly, the sensitivity to is given by: 


C f 

3K^ 

5^ ^ ~ 025 

K. c f 

a a 

Because of the very small floating tendency of the flaps, the 
aircraft is virtually insensitive to uncertainties in the flap 
hinge moment derivative with respect to changes in angle of 

attack. Thus, a 100% error in c£ only produces a 2.5 % uncer- 

. ot 

tainty in K^. 

The same cannot be said, however, of the effect of uncer- 
tainties in the determination of the aerodynamic stiffness of the 
flap. Since the auxiliary wings do not contribute to the overall 
aerodynamic stiffness of the system, K is inversely proportional 


to C 


Hi 


Thus, a 1% underestimation of C„ results in the same 

“5 


percentage overalleviation. Although the uncertainty in the 
theoretical estimate of this derivative can be reduced to some 
extent by experimental measurement of the appropriate hinge mo- 
ment coefficient, it must be known within 10% of its exact value 
if is to remain within the acceptable range of 0.9 to 1.1. 

However, the uncertainties in C„ theoretically can be com- 

“6 

pensated by varying the auxiliary wing size L or the gearing 
ratio r. If the auxiliary wing aspect ratio is kept constant. 


- 9K K 

— — ^ = 3(1- — ^)=!3 

9L ' K ^ 

a a 


(32) 


If only the span is changed, the dependence of the lift curve 
slope on the aspect ratio causes the relative sensitivity to be 
slightly in excess of 2. 


L 

K 


a 



9L Cf" „f 

n L K 

n 3 . 


(33) 
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It would appear that one could make use of the large sensitivities 
of K to variations in auxiliary wing size to overcome uncertain- 

cl 

ties in C„ . However, it is rather impractical to trim the aux- 
^ 6 

iliary wing dimensions and, fortunately, adjustment of the gear- 
ing ratio offers a more realistic compensation technique, and is 
fairly easy to implement. 



(34) 


Thus, the absolute sensitivities of K_ to the gearing ratio and 

a 

to C„ are the same. 

In order to trim the aircraft to the full vertical allevia- 
tion condition, it first should be flown with a small enough gear- 
ing ratio so that it is effectively in an underalleviated con- 
dition (with the pushrod attached fairly close to the flap hinge) 
The gearing ratio then should be increased progressively by in- 
creasing the flap moment arm until the time to halve the pitching 
response of the aircraft to a step of elevator increases to about 
0.9 second, corresponding to that for the flap-configured Cessna 
172. This procedure accommodates small uncertainties in C 

m ^ 

and 3e/36, which have a marked effect on the short -period natural 
frequency but relatively little effect on the time to halve ampli- 
tude. As a final check on the value of K , the phugoid oscilla- 
tion should practically disappear from this response, the indi- 
cated airspeed should remain constant and the altimeter should 
show no short term perturbation. This entire test program, of 
course, should be undertaken without having the elevator coupled 
to the auxiliary wing. Finally, the incidence of the auxiliary 
wing should be adjusted so that the flap remains undeflected in 
the absence of gusts. 

The following analysis investigates the effect of flap 
pitching moment and downwash derivatives on the short period dynam 
ic stability of the vertically alleviated airplane. For K^=0 
and K_=l the short -period natural frequency and damping ratio are 

a 
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determined by; 
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The augmented pitch stiffness, C takes into account the modi- 
fied lift distribution of the wing and the change in the down- 
wash at the tail resulting from flap deflection. 


where 
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The short -period roots can be traced in the complex frequency 

d £ 

plane as a function of and , yielding two sets of curves. 

""Sw 

3 £ 

As shown in Figure 21 , curves of constant flap downwash, are 
simply straight lines parallel to the imaginary axis while those 
of constant flap-wing pitching moment are slightly curved. The 
intersection of two members of different families yields a pair 
of complex conjugate roots representing the oscillatory mode 
corresponding to a particular set of flap characteristics. If 
the two curves do not intersect, the short period mode is de- 
generated into a pair of real roots which are symmetrical about 
the corresponding ^ line since remains the arithmetic mean 

of the roots. It should be noted tEat the amplitude of the pitch 
ing response to vertical gusts of high frequency is equal to 

twice the distance between the roots (or the mean of the roots 

D£ 

if the mode is degenerated) and the line = +1 whrch 
represents optimal vertical gust alleviation. This distance 
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rad/sec 


Figure 21 Effects of Flap Downwash and Pitching Moment 
on the Short-Period Roots of the Vertical- 
Gust Alleviated Aircraft 
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is simply: 

a = |(i) (39) 

(jO-9-oo 

Several comments are in order. By somehow decreasing the magni- 
tude of the adverse nose-down pitching moment of the flaps, sub- 
stantial increases in the natural frequency of the short-period 
mode can be achieved without changing the high frequency pitching 
response to gusts or time to halve amplitude of the damping ex- 
ponentials. For example, reducing C to zero shifts the roots 

'^w 

from point (a) to point (b) in Figure 21 and increases the natur- 
al frequency by a factor of four. This reduction in flap pitching 
moment in principle could be obtained by proportional coupling 
of the elevator (or a portion of the elevator) to the flaps so 
that an upward gust would produce up flap and down elevator de- 
flections. Then if the flap downwash derivative could be reduced 
until the damping of the mode increased to the desired level, 
substantial reductions in pitching response could be achieved 
while also improving the stability characteristics of the alle- 
viated aircraft. For C =0 and ^ =.5 corresponding to l-f- =-.022 

m^ Da ^ do 

°w 

(point (c) ) , the pitching and load factor responses to high fre- 
quency gusts can be reduced by a factor of 3 compared v;ith point 
(a) while the natural frequency is approximately doubled and the 
damping ratio is increased to .5, which is the same as the basic 
aircraft. 

It should be pointed out that flaps with < 0, that is 
creating upwash on the stabilized when lowered, produce an 
attenuation in the pitching response relative to the basic air- 
craft (^ = .37, which is to the right of point (c) in Figure 21). 

If the flexibility of altering the flap geometry is not 
available to the designer and elevator coupling is too complex, 
he is left with the option to couple the ailerons to the flaps to 
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reduce the magnitude of the downwash while maintaining a stable 
aircraft. The results of this process are given in Figure (22) 
which shows the locus of the short period roots of an alleviated 
aircraft in which the gearing between the flaps and symmetric 
mode ailerons is varied. As can be seen, reducing the downwash 
derivative 3 e/3 6 by positive coupling (flaps and ailerons in 
phase) reduces the pitching response slightly but quickly leads 
to degenerate and unstable modes. Because of the strong adverse 
pitching moment of the flap, the range of permissible aileron 
coupling is sever ly limited and offers little advantage. 

However, test data for a model of the Cessna 172 generated 
too late in the program to be used in any but the sensitivity 
analysis are applied to Figure 22 to provide a second 
locus as a function of aileron coupling, which provokes consider- 
ably different conclusions. The very low tail-off pitching mo- 
ment derivative with flap deflection and a slightly smaller value 
for 3e/36 yields a curve which in essence is shifted upwards by 
about 4 1/2 units and to the right by one quarter unit. Because 
of the resulting higher frequencies, the sensitivity of the roots 
to aileron coupling is reduced substantially. Most importantly, 
aileron coupling not only becomes beneficial, but the ailerons 
should be considered the principal control surfaces and moderate 
flap coupling is helpful in avoiding the real root with long time 
constant introduced by using the ailerons and to yield short 
period roots that satisfy handling criteria. 

With the use of ailerons, the downwash affords a more favor- 
able high frequency pitching response. If the aileron-flap 
coupling ratio, G' , is 5, there is no change in this portion of 
the response when the vertical gust alleviation system is imple- 
mented. Even with no coupling, the response is only doubled, 
compared to the tripled value for the analytical model. 

It should be noted that test data are not available for the 
aileron parameters. These were calculated by methods presented 
in Reference 6, as were the flap parameters. Since the tail-off 
flap pitching moment derivative with deflection was twice 
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Figure 22 Effect of Aileron Coupling on Short 
Period Roots 
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the magnitude of the measured value, a similar relationship was 
assumed for the aileron derivative. The test value of C 

w 

given in Figure 22 was derived by scaling the value computed at 
the quarter m.a.c. (-.235) by the ratio of tested and computed 
flap derivatives at this point (.25/. 485) and then referencing 
the derivative to the 33 percent center of gravity location by 
making use of and the eight percent m.a.c. moment arm. 

^a 

In conclusion, introduction of the test data emphasizes^ the 
need for accurate measurement of the flap/aileron characteristics. 
Short-period roots are quite sensitive to uncertainties in these 
parameters, particularly to C . If the full scale aircraft 

lll^ 

w 

reflects the model data, the short-period mode would be lightly 
damped if flaps are used to alleviate gusts. If ailerons are 
employed, the roots are real and stable. Use of the ailerons 
with modest coupling of the flaps would yield roots preferable 
with regard to handling qualities. 
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VIII . CONCLUSIONS AND RECOMMENDATIONS 


A passive aeromechanical gust alleviation system is developed 
(specifically for a Cessna 172) using existing lift and control 
surface regions. 

Conditions for optimal gust alleviation are developed and 
shown to be incompatible with normal stability requirements. A 
stability analysis in the form of root loci is conducted to iden- 
tify acceptable flap characteristics. 

The efficiency of the gust alleviation system is examined 
in terms of load factor and pitching responses to gusts over the 
broad frequency spectrum. The Von Karman model of atmospheric 
turbulence is then applied in order to investigate the relative 
importance of these responses at the short period and phugoid 
frequencies. As predicted, the considerable attenuation in normal 
acceleration is accompanied by a substantial increase in pitch- 
ing response. The effects of the alleviation system on the 
maneuverability of the aircraft are then analyzed and the concept 
of coupling the pilot's stick to the flap system is developed to 
explore the possibility of recovering some of the controllability 
normally associated with elevator movement. Finally a sensitivity 
analysis is conducted to determine the designer's ability to com- 
pensate for uncertainties in knowledge of flap/aileron hinge 
moment characteristics. 

The analysis has shown that passive gust alleviation is 
practical. Gust loads can be absorbed well at frequencies be- 
tween the phugoid and short-period mode where the important gust 
energy is concentrated. The major task is to determine the flap/ 
aileron characteristics on which to base the auxiliary wing design 
In order to minimize the size and additional weight of the sensors 
the lift control surfaces should be properly balanced about their 
hinge line in order to provide the largest possible aerodynamic 
floating tendency. Also a fast flap-auxiliary wing response is 
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necessary to achieve good alleviation at higher frequencies 
(near the short»period and above) . Thus the system should be 
made as light as structurally feasible while also maintaining 
the damping ratio within acceptable bounds. For a given 
alleviation constant (but variable gearing ratio) the inertia , 
and aerodynamic damping (referred to the flap hinge line) con- 
tributed by geometrically similar auxiliary wings are inversely 
proportional to their size and size squared respectively. 

Thus the sensors should be made as large as possible subject to 
limitations of weight, drag and linkage flexibility. 

The pitching moment and downwash characteristics of the lift 
control surfaces are also shown to have a large effect on the 
aircraft's response to gusts above the short-period frequency. 

The pitching motion of the aircraft can be aggravated by the 
alleviation system if the nose down pitching moment of the flaps 
due to increased camber of the wing is excessive. Alleviation 
of this motion as well as improved load factor attenuation can 
be achieved if the adverse flap pitching moment is small enough 
to allow substantial reduction in the downwash produced by lower- 
ing the trailing edge surfaces without jeopardizing the aircraft 
stability. This reduction in downwash can be achieved in practice 
by coupling symmetric mode ailerons to the flaps. 

Maneuverability is virtually eliminated by the alleviation 
system but can be fully recovered by proper auxiliary wing-ele- 
vator coupling and even enhanced by the effective introduction 
of direct lift control. 

Finally, residual uncertainties in flap hinge moment parameters 
can be accomodated by adjustments in the flap-auxiliary wing gear- 
ing. A safe and practical procedure for implementing the alle- 
viation system is presented. 
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RECOMMENDATI ONS 


In view of the importance of flap-auxiliary wing parameters 
on the performance of the gust alleviation system a general 
parameter optimization program could further enhance the design 
process. Also the low speed handling characteristics of the air- 
craft at take-off and landing need to be examined in order to de- 
termine whether or not the alleviation system should be disengaged 
and if necessary to design a cockpit activated decoupling or 
locking device. 
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APPENDIX A 


Cessna 172 Par^ters 

The Cessna 172 parameters used in this study are presented 
in the table below. They were derived by the procedures given 
in the Air Force DATCOM, except for the rotary derivatives which 
were taken from NASA CR-1975 in which Cessna 182 parameters were 
developed. These aircraft are sufficiently alike to warrant the 
use of these data. Very near the end of the program, NASA fur- 
nished data based on wind tunnel tests of a six-foot span model 
of the Cessna 172. This provided measurements of flap character- 
istics which are difficult to estimate, and confirmed a number 
of the computed stability derivatives. 


TABLE A-1 

Cessna 172 Parameters 



Program 

NASA 


0.416 


Cl 

5.50 

5.50 

a 



Cl 

5.07 

5.03 

“w 



Cl 

0.66 

0 . 62 

“t 



3e 

3a 

0.37 

0.25 


1.49 



3.88 


*c„ 

-0.83 


*C 

m* 

a 

-4.36 


*c.« 

-11.40 



77 



Line 



8 





TABLE A-1 cont 



Progreun 

-1.26 

0.25 

-0.116 

99.9 

287 (1581 Kg-m^) 
59.13 m/sec 

1000 kg 
2.30 m 


Flap Characteristics 



flap 

aileron 

flap-aileron 

NASA Test 
(flap) 

^^6 

w 

1.03 

.75 

1.78 

1.10 


.598 

-.140 

.458 

- 

Ift C 

- .403 

-.175 

- .578 

-.168 

*C 

mj 

-1.75 

.408 

-1.34 


de 

VS' 

.150 

-.035 

.115 

. 15 


.034 

.034 

.034 

- 


-.46 

- 

— 

-.40 


-.068 

- 

- 

- 


* Values shown for CG at 33% MAC, static margin 15% MAC 
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The parameters 


m. 


9 e 

and ^ have approximately linear 


w w 

relationships among the configurations; that is, flap values 
plus aileron values equal flap-aileron values. (Ref. 6, Sec. 4. 4 
and 6.1) . The flap downwash lag derivatives are assumed due 
only to tail loads and are found from the formulations given in 
Appendix E. As a comparison with the DATCOM, the flap induced 
downwash at the tail also was computed from NACA TR 648. The 
derivative 9 e/9 6 was found to be .25 for the inboard flap and 
.178 for full span flaps. Use of the DATCOM-derived figures 
assures conservative stability estimates. 

The flaps and auxiliary wing physical parameters are given 
below, and are followed by the computation of the system natural 
frequency. The mass balance for the flap has been neglected in 
computing the flap moment of inertia. However, its contribution 
to the moment of inertia of the system is negligible. 


MOMENT OF INERTIA OF FLAP 
ABOUT ITS HINGE LINE 


1) Spar 

Weight = .33 kg 

M.I. = .0037 kg-m^ 


2 ) Ribs 

Vfeight = . 31 kg 

M.I. = .0079 kg-m^ 


3) Leading Cover 

Weight = 1.27 kg 

M.I. = .0120 kg-m^ 
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II 

4) Trailing Cover 

Weight = 1.65 kg 

M.I. = .120 kg-m^ 

Weight of Flap = 3.57 kg 

,M.I. of FLap = -.143 kg-m^ 


MOMENT OF INERTIA OF THE AUXELIARY WENG 
ABOUT ITS HINGE LINE 


1) Spar 


Weight = 

M.I . 


. 12 kg 
.053 kg-m^ 


2) Ribs 

Weight = . 11 kg 

M.I. = .056 kg-m^ 


3) Leading Cover 

Weight = . 55 kg 

M.I . = .232 kg-m^ 


4) Trailing Cover 

Weight = 1.21 kg 

M.I. = .51 kg-m^ 


5) Linkage 

Weight = .36 kg 

M.I. = .034 kg-m^ 


6) Counterbalances 

Weight = 7.10 kg 

M.I. = .25 kg-m^ 

Weight of the Auxiliary Wing = 9.45 kg 

Moment of Inertia about Hinge Line = 1.14 kg-m 
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The moment of inertia of th^ flap-auxiliary wing system, 
about the . flap hinge line depends upon the gearing ratio. 


If,aw= 


(A-1) 


The damping term in the hinge moment equation is primarily a 
function of the "roll-damping" developed by the auxiliry wing. 
Simple strip theory can be used tP estimate this coefficient.. 
If X represents the spanwise ordinate along this wing, 


H 


aw 


I 

-/ 


aw ^ 


=- ^(Sb2)^^ c^ 

ot 


(A-2) 


(A- 3) 


where Ct is the section lift curve slope of the auxiliary wing, 
■“a 

The moment and angle transmitted to the flap hinge are: 


6 


aw 


r 6, 


H = H_ + TH ' 
f a\ 


Then, 


'H; 


1 3H 

q(Sc).p 5 -c 


= C 


H; 


o o (Sb) ,aw 
2y2 ' aw b 

3 (Sc)^ c I 


or, 


C„ --.6 - 6.2 = -6.8 
^6 


(A-4) 
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A more accurate result is obtained particularly in the case of 
low aspect ratio vanes if the expression for the roll d^unping 
of a wing is used instead. 




(A-5) 


where Cg^ is the roll derivative based on twice the auxiliary 
wing spaH. 

This method yields a much smaller value for the dcunping contri- 
buted by the vane, 

C„. = -.6 - 1.9 = -2.5 
“6 

From the foregoing analysis and consideration of the flap de- 
flection stiffness, the system characteristic equation can be 
written. 

The natural frequency of this system is given by: 

= 0.166 (13.2 rad/sec) (A-7) 



where : 


Inertia , 1.28 kg-m^ 

c^ 2 

p(Sc)^.'g- .0761 kg-m 


16.8 
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APPENDIX B 


Analog Simulation of Aircraft Longitudinal Dynamics 

Simulation of basic, gust-alleviated and control-recovered 
aircraft longitudinal dynamics was conducted on a small, 10 volt 
analog computer in order to investigate the response in real time 
of the modified aircraft to gust and control inputs. The flap 
avixiliary-wing system dynamics also were included in the simula- 
tion in order to identify any adverse effect that might occur as 
a result of interference with the short-period mode of the air- 
craft. The high inertia of the initial design resulted in only 
50% load factor attenuation at the basic aircraft short period. 
However, by increasing the natural frequency of the system and 
reducing the aerodynamic damping correspondingly, the idealized 
levels of attenuation calculated in the analytical treatment 
could be reached. 

An attempt was made to simulate atmospheric turbulence by 
feeding the output of a"white" noise generator through a low pass 
filter in order to reproduce the Dryden spectriam. Unfortunately 
lack of a noise generator with good low freqency characteristics 
made it unsuccessful. Discrete gusts were simulated using a well 
damped second-order system. Gust shape and duration were adjusted 
by controlling the damping and natural frequency of the double in- 
tegrator feedback system. 

A schematic diagram identifying the scaling factors (s^, Sg, 
etc. ) and corresponding analog quantities in brackets is shown in 
Figure B-1. Elevator deflection was scaled at 100 volts/radian. 

All other angular quantities were scaled at lOv/rad. The nondimen- 
sional velocity perturbation was scaled at lOv/unit. t=^-. 0125 sec. 
designates the non-dimensionalizing time sealing factor. 
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APPENDIX C 


On the Mechanization of the Gust Alleviation Device 


The degree of vertical gust alleviation as measured by the 


general alleviation constant is determined by: 


'H 


^a = C 


a 


H, 


(C-1) 


‘6 "a 

The free floating aerodynamically unbalanced flap does not con- 
tribute much to the improvement of riding qualities in turbulent 
air since it only yields a very small value of whidh is due 
to its relative insensitivity to angle of attack gusts, 

, C„ C, 


- / H L.\ 

^ -.025 


(C-2) 

c ' 

0 a 

As demonstrated previously, full vertical gust alleviation re- 
quires to be unity. 

Means of increasing the flap response to angle of attack 
(without the use of servo-mechanisms) essentially fall into 
two categories : 

Increasing the flap hinge-moment sensitivity to angle of 
attack C„ (by a factor of 40) by coupling a pair of auxili- 
ary wings to the lift control surfaces in order to provide 
enough power to drive them. This type of device was adopted for 
the present system. 

Reducing the coefficient (by the same factor) in order to 

increase the aerodynamic balance of the flaps to the required 
level. This could be done in two ways; 

1. Mechanically, by providing a destabilizing spring 
moment about the flap hinge. 

2. Aerodynamically, by installing a servo- tab at the trail- 
ing edge of the flaps. 

The spring mechanism would appear as shown below: 
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The spring offset angle 3 is introduced to provide horizontal 
gust alleviation capability as measured by K^. The hinge 
moment provided by the spring about the flap hinge is: 

H = Fr sin (6+3) (C-3) 


The tension F in the spring can be made constant by providing a 
large enough equilibrium extension. For zero flap deflection. 


and 


H = Fr sin 3 
eq 

H- = Fr cos 3 


(C-4) 

(C-5) 


The resulting hinge moment coefficients for a particular flight 


condition become: 

Cl 

n + Fr cos 3 ^ 1 6 

V «6 

Cfl “ 

_ Fr sin 3 _ p 6 

q(Sc). C, 

r f Lg 


or. 


tan 3 = 



( ^ 


1 


) 


'H 


a 


(C-6) 

(C-7) 


(C-8) 


Thus if only vertical gust alleviation is to be implemented K =0 
and the spring offset angle 3=0. However, even if full vertical 
and horizontal gust alleviation are introduced, 3 is only .66 de- 
gree. The magnitude of the required spring loading is determin- 
ed from equation (C-6), 


Fr cos 3 = -q(Sc)^ 




= 160 ft lb. (C-9) 
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Two comments are in order: 

- At the most r is of the order of 1 ft with the consequence 
that a rather large transverse force is applied to the flap 
hinge and may cause structural problems. 

- The spring tension must be trimmed to the correct value for 
the equilibrium dynamic pressure even if only vertical gust 
alleviation is implemented. 

These problems are eliminated if a servo-tab is used to improve 
the aerodynamic balance of the flaps. 



The tab is geared to the flap 
flap is lowered thus reducing 

C = (1 + 

“6 

in order to stay in the range of linear aerodynamics, the gear- 
ing ratio should be reasonably close to unity: 


so as to deflect upward when the 
the magnitude of C 


H, 


to : 


96 ^ 


'H, 


'H, 


-) C 




(C-10) 



(C-11) 


The required tab moment Cjj 


1 + 



is determined by the condition: 



(C-12) 
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Thus, for full vertical gust alleviation, the tab should 
be designed so that 



If horizontal gust alleviation is to be included, the flaps must 
be spring loaded as previously considered or the tab can be 
spring loaded. 

The elimination of the auxiliary wing results in a "clean- 
er" and probably less expensive alleviation device. Also the 
inertia of the system is greatly reduced when the auxiliary wings 
are removed (by a factor of 9) . However the reduction in the 
restoring moment is even greater (1/40) so that the natural 

frequency of the system drops to: 



: 6 . 3 rad/s 


which is approximately half the present value. This value may 
interfere with the alleviated aircraft short period and intro- 
duce instability. However, further analysis shows that a com- 
bination of a spring and an auxiliary wing can raise the system 
natural frequency above that for a system with an auxiliary wing 
alone. 

In the presence of the auxiliary wing 



,f Ij, ^^^^aw -aw 

"c. ■ 2 <So)f X 


(C-13) 


where T, 


35 

T6 


aw 

f 


In the presence of a spring or tab: 

- hJ + Pr 


(C-14) 
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For full vertical gust alleviation, K_=l , 

cl 




a 


a 


'H, 


(C-15D. 


‘6 "6 

For geometrically similar auxiliary wings: 


I - l""®* 
aw aw 


(C-16) 


where the reference auxiliary-wing characteristics are those 
already calculated: NACA 0009, AR=2, bJ^®f=1.13m , lj®^= 1 . 14kg-m^ 
rref^i, (l£=.143 kg-m^) 

After some manipulation, equations (C-13) , (C-14) ,and (C-15) yield the 

H 

6 Fr 

following relation between parameters r, L, and X= .= — =1+ = — 

n. jf il 


- X-K 

TL^ = - a 


1-K 


"f 

(C-17) 


The natural frequency of the alleviation device is given by: 


0,2 = - 

n 




(C-18) 


from (16) , (17) and (18) , 


/ 0) V 2 

(■c) • 


.L 

r 

a 


/X-K^ / I. 


where w is the free floating flap natviral frequency. 


n 
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For small values of K^: 

a 



X 


1+X 


iref 

aw 


Then 


0 ) 


n 


w 




The ratio 



is just a function of 

has been plotted against the spring dependent 


parameter X yielding a family of curves corresponding to differ- 
ent auxiliary wing sizes. For a given auxiliary wing size there 
is a value of the spring moment corresponding to a maximum 
natural frequency of the system. For all practical auxiliary 
wing sizes this optimum moment is destabilizing. The limiting 
size corresponds to L = 8.4 for which no improvement 

in natural frequency can be obtained through the use of a 
spring; but this case is rather academic since the auxiliary 
wing span would then be approximately twice that of the wing I 
For a given spring moment (which may be zero) the larger 
the auxiliary wing size, the larger the natural frequency of 
the system. For the present auxiliary wing, b^'^=3.7ft (L=l) , 

Xmax~*845 corresponds to Fr - 107 ft lb and / ‘^n \ ^ = i X 

lo) / 2 max 

nf/ 

or; w =16.8 rad/s since to =40.4 rad/s 
n nj; 

max f 


For 


X = L = 1 
No spring 


For 



( No aux.wing 



.326 


= 13.2 rad/s 


0 )^ = 6.3 rad/s 
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Figure C-1 Variation of Flap-Auxiliary Wing Natur 
- of Restoring Moment 


positive restoring moment 
spring or tab 

Frequency as a Function 




APPENDIX D 


Simplified Frequency Analysis 

I. Analysis of Frequency Responses to Gusts 
The Short Period Approximation 

The velocity perturbation is assvuned to be negligible, with 
the result that the equations of motion simplify to. 


pi, + 

s(C^^ + 2y) 

s (C. 

a 

a 


C* + 

s C ’ 

s(C 

m 

a 

m- 

a 


mq y 


a 


e 


a a q 

C* + s(C* “ C ) 
m m • m 

a a q 


a =0 

«g 


(D-l) 

If it is assumed that the tail is the only contributor to the 
pitching rate and angle of attack rate stability derivatives, the 
relationships given by Phillips and Kraft^ yield the following 
identity. 


Ct C - C = 0 

m • m„ L* 

q a q a 


(D-2) 


The angle of attack transfer function then simplifies to; 


a 

a_ 


C -C 2 — — — 2 

s^ + 2? 0 ) s + oj ^ 


a 


sp sp 


(D-3) 


where 


(0 2 
n 


sp 


- 2 
“n 


c; C + (2y-C^ )C 

"^q ^q 

I (2y+>C'L.) 

^ a 

Cl C + (2y-C_ )C' 

L m m 

a q q a 

q g 


Cl I - 2y(C +C'.) 

L y m m* 

2CW = « 9 «- 

%p I (2y+C^J 

g 

Ci I + 2y(C -C ) 
L y m m* 

25 S 3-L a s_ 

”sp I (Cj, -Cl.) 

q g 
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similarly for the pitch angle. 


2y.(C - C )s + (2y+C_. )C -C <C1 

^ m*- ' m L m • m ■ L 

0 _ a g q a q a 


a. 


I (2y+C' . ) (s^+2c co_s + u 

y ■ . ^sp "sp 


a 


The flight path angle response then becomes. 


Y _ 9 -.g _ 


-a 


C -C ' 

' L l: 

q g 

2M-fC£. 

a 


2 - - -2 

s +2^ 0) s + 0) 

n n 

2 2 

s +2? tjo s +' 0) 
n n 

sp sp 


where 


to 


n 


2 (c c; - c' ) 

m L L m 

q a q a 

Iy<=L 

q a 




2? w = - 


a 


n 


c -c * 

L* 

q a 


2- ri -A 

(D-4) 


(D-5) 


2 2 

It should be noted that terms of the form (s +2Cco s+to ) dp not 

n n ' ^ 

necessarily represent oscillatory modes of natural frequency 

and damping ratio 5 , but are general second-order terms that may 

very well be nonoscillatory . Such is the case with the zeros of 

the basic aircraft's flight path angle transfer function where 
2 

60 ^ , is negative and the term has two first-order roots; one is 
large and positive at approximately 

s 

i 

and the other small and 

s 


a 


L L* 
q a 

negative 

2C 

m 
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The Phugoid Approximation 


The phugoid approximation is based on the assumption that 
the total angle of attack perturbation (a+a^) is negligibly small 
and that the pitching moment equation is identically satisfied. 
With these assumptions the equations of motion reduce to: 


c* - 2ys 

-^L 


vT 


0 

*u 






2C- 


• 


+ 

) 

2Cj,U-K^)- 

L ^6 ® 

s(C - 2y) 


e 


-s(C +2y) 

q 


— - 


a _ 


a =0 

g 


(D-6) 


The velocity perturbation transfer function is then: 


s 


u 

“g . ^ 0.2 

P P 


(D-7) 


Similarly for pitch angle. 


s-C' /2y 
*u 


ot. 


p p 


• s 


(D-8) 


The flight path angle response thus becomes. 


s^+2^ 0) 0) ^ 

_ ^p n„ Ct n^ 

Q L ^ p L p 


X = (_H_ +1). a 

“g “g V 


2 2 
s +2 5_ti>_ s + 

P 

P P 


(D-9) 


where : 


u. 


2 ~ 
n_ 


2y^ 


C 

X 




u 


2y 
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Ij 

For the basic airplane, (o^= yT?' phugoid frequency 

range the flight path angle reduces to: 





n 


2 

P 


* “n„ 
P P 


(D-10) 


For the airplane equipped with full horizontal gust alleviation 

capability, K^*lf =0, and 

P 



(D-11) 


The simplified frequency response analysis permits the con- 
struction of approximate Bode diagrams of various quantities. 
Attention is concentrated on the load factor and pitching rate 
which are of particular interest with regards to riding qualities 
of the aircraft in turbulent atmospheric conditions. 

The non-dimensional rates of change of flight path angle 
and pitch attitude are related to the corresponding dimensional 
load factor and angular pitching velocity response; 

SI (g-sec/m) 

a V gc 'a„ ' 

g g 

^ . I (^ ) (raa/m) 

g g 

since the effectiveness of the gust alleviation system depends 
essentially on the level of attenuation of the load factor re- 
sponse in the short-period mode frequency band, it is important 
to be able to predict the response of both the basic and the alle- 
viated airplane in that range of frequencies. 

The basic airplane is seen to have a flat load factor re- 
sponse between the short-period frequency u) and the upper break 

sp 


(D-12) 

(D-13) 
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I 


frequency at 



A simple extrapolation of the high frequency asymptotic 
response down to the upper break frequency enables us to determine 
its level. 


JL cs - 


a 


g 



or 


An 

2V . 

a 

a V 

g 

~ 0 

gc 

2y 


.223 

(g-sec/m) 


(D-14) 


The basic aircraft's pitching response also levels off beyond the 
short -period frequency; 


- 2 . 


a 


g 



or 



2 

c -c 

m* m 

a 

ca_ 

c 

I 

g 


y 


.031 

rad/m 


(D-15) 


In a similar way the high frequency response of the alleviated air- 
craft can be approximated. The load factor's asymptotic behavior 
is characterized by: 


JL == 
a 

g 


Ct -c^ 

L L* 

g a 

2y 


or 


An 


. q a 

a V 


2y 

g 



- =^3.02 10“^u 
^ (g-sec/m) 


(D-16) 


where w is the gust frequency in radians per second and is relat- 
ed to the nondimens ional Laplace variable s through s=io)^ with 
i=/^. 

The pitching motion response behaves similarly to that for the 
basic aircraft beyond the short period frequency. 



= .082 
rad/m 


(D-17) 
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A similar analysis has been made for frequencies below the phugoid 
Resulting expressions are shown on the Bode plots in Chapter IV 
along with those presented above. 


II. Responses to Elevator Deflection 

Under the assumption of instantaneous flap deflection, the 
response of the airplane in horizontal cruising flight to elevator 
deflection is found by solving the following matrix equation: 


fc' -2ys 

^U 

c 


-^L 


U 


C 







I4'=I./==L5> 

ot 

+ s(2y+C^J 
a 

s (C -2y) 


a' 

+ 


(c +sc ) 
[su 

C* 

m 

a 

+ sc* . 

m* 

a 

s(C -Is) 

iHq y _ 

■ 

0 


C 

m. 


(D-18) 

The preceeding equation was used directly for machine computation 
of frequency response functions which basically require only rou- 
tine operations on matrices with complex coefficients. The result 
ing velocity, angle of attack, pitch and flight path angle pertur- 
bation amplitude Bode diagrams are presented in Figures 16 - 19. 

However, for analysis one needs literal expressions for the 
various transfer functions which is impractical for the complete 
system especially if the inverse Laplace transforms are required 
to study the system in the time domain. The exact solutions, in- 
cluding the effect of flap dynamics were obtained using real time 
integration of the original differential equations on an analog 
computer . 

For analytical work associated with control system design, 
approximate forms of the transfer functions are useful and are 
readily obtained by extension of the previously developed short 
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perip4 and phugoid approximations to the case of the elevator in- 
puts. 

The Short Period Approximation 

The short period mode of the longitudinal motion is essen- 
tially a constant: speed oscillation, in which normal and pitch- 
ing motions occur. The characteristics of the mode, that is the 
frequency and damping, together with the control system dynamics 
and sensitivity to gusts, define the short term response to con- 
trol inputs and to gusts, and therefore can have a large influence 
on a pilot's views of the handling qualities of the aircraft. 

It is assumed that u^O and that the speed equation of motion 
is identically satisfied. 

The system of equations thus reduces to: 


C' +s(2y+C' ) s(C -2y) 
a ■^a "^q 


C' +sC 

01 a 


s(C -Is) 

q ^ _ 


/a\ 


6 =0 
e 


(D-19) 


m. 


Under the assumption that led to the determination of stability 
derivatives involving tail lag effects. 



The angle of attack transfer function reduces to: 





°e 

2 

s + 2^01 s+o) 
"^sp 


6 

2 

n 


e 


sp 


(D-20) 


(D-21) 
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while the pitching angle reduces to: 



s+(C; C -C’ C_ )/2yC 

ot 0 ^ a 0 ^ 0 

e e e 


s (s + 2 ^ 0 ) s+o) ) 


n 


sp 


n 


sp 


(D-22) 


The flight path angle then follows: 


JL = 

^e 


2V 


s%(c; c -c' c )/i 

e e e 

2 5 

s (s +2^0) s+(i) ) 

n n. 


(D-23) 


sp 


sp 


The short period mode natural frequency and damping ratio were 
determined previously and are given by: 


b). 


n. 


sp 


C' C +2yC 
L m m 

g q g 

2yi„ 


2cu) 


sp 


C’ I -2y(C +C ) 

^c. y % % 

2UI„ 


In view of the fact that the zero appearing in the angle of 
attack transfer function. 


2yC 


m. 


s.= 


g I. C 


y Li 


-162 rad/s 


is very large compared to the short period natural frequency, the 
angle of attack perturbation is well approximated in the frequency 
band of interest by the following transfer function: 




s ^+250) 


n 


S+(i) 


sp 


sp 


(D-24) 


Between the phugoid and short period frequencies the angle of 
attack response is flat, the level of the response being strongly 
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dependent on the degree of vertical gust alleviation as determin- 
ed by K . 


VO (Basic) ^ - c, C *ivC = 
e Xj m m 

a q a 

C 

"*6 

K,=l (Alleviated) #- = - 7 ^ 1 —^ = -8.67 

^ ^e 

a 


(D-25) 


(D-26) 


The pitching response in the short period frequency band is domin- 
ated by the presence of the finite zero, 
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the position of which is strongly dependent on K . Thus, for the 
basic aircraft, K_=0 
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The initial attitude response to a step input in elevator is 
therefore well represented by a ramp corresponding to the steady 
pitching rate. 
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On the other hand, for the alleviated aircraft, K =1, the zero is 

a 

negligible compared to the short period natural frequency. 


C C- 
m L. 
a 0, 


Sfl = 


2y C, 


m. 


0.02 rad/s 


The pole at the origin is therefore effectively cancelled and the 
pitch angle reduces to: 
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The initial attitude response is that of a simple second order 

system of natural frequency co :1.82 rad/s and damping ratio 5 =.41 

n 

with a steady- state gain 


-^t-^=-8.67 (D-30) 

® s ->0 "‘a 


It should be noted that Within the limits of the approximations 
made, the aircraft equipped with full vertical gust alleviation 
has identical pitch and angle of attack transfer functions in the 
short period frequency band so that the flight path angle and 
consequently the altitude are not affected by elevator deflections 
at these frequencies. 

If we return to the original expression of the flight path 
angle, the basic aircraft response is given by; 
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whereas when full vertical gust alleviation is implemented, 
the response is strongly attenuated. 
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Thus, in the short period frequency band, the alleviated aircraft 
effectively "climbs on it's tail" since 



(D-33) 


except at the modified short period natural frequency where an 
undamped pair of oscillatory zeros cause the vertical response 
to drop sharply to zero. Note the reversal in control. 


The Phugoid Approximation 

The phugoid is a lightly damped low frequency oscillation in 
which changes of airspeed and height occur at almost constant 
incidence. When the elevator angle is varying however, the angle 
of attack perturbation can no longer be assumed small. The equiva 
lent assumption in this case is that pitch equilibri\am be maintain 
ed. 
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At low elevator frequencies, the above expression can be replaced 
by: 
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The spring loading of the flaps is designed to attenuate the load 
factor response to variations in dynamic pressure resulting prim- 
arily from horizontal gusts. It also will make the flaps sensi- 
tive to changes in airspeed induced by control deflections or angle 
of attack gusts exciting the phugoid mode and thus drastically 
modifies the basic aircraft's low frequency characteristics. 

More specifically, implementation of the device virtually 
eliminates the load factor response to horizontal gusts over the 
broad frequency spectrum, up to 10 rad/s. In the cruising con- 
figuration, the sensitivity to horizontal gusts is relatively unim- 
portant. However, when landing at low speed and high lift coeffi- 
cient, horizontal gusts have a greater influence on the normal 
acceleration. Furthermore, the phugoid frequency is higher at this 
speed and may be more difficult for the pilot to control. Sub- 
stantial attenuation of the normal acceleration in low frequency 
vertical gusts and in particular suppression of the resonant 
peak at the phugoid mode is another benefit of spring loaded flaps 
although of minor significance since, in general, the pilot 
should have little difficulty in controlling that mode. There is 
however a fundamental deficiency in the stability characteristics 
of an aircraft equipped with a gust alleviation system designed to 
satisfy the constant static lift condition which is immediately 
apparent from the lift equation and results in marginal stability 
of the airplane. This makes controllability unacceptable without 
adequate corrective measures. 

With these facts in mind we will first concentrate on the 
phugoid approximation as applied to the basic aircraft and to the 
alleviated aircraft without a spring restrained flap system. 

With K^=0 and at low elevator frequencies, the pitch equation 
requires that the angle of attack be proportional to the elevator 
deflection: 
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So that- the steady state angle of attack is approximately multi- 
plied by a factor of 6. Substitution in the remaining lift and 
drag equations yields the pitch angle transfer function. 
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while the velocity perturbation is given by 
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Expanding and simplifying the numerator determinant yields the 
basic aircraft velocity transfer function: 
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The derivative term can be neglected at low elevator frequencies 
because the zero is large compared to the phugoid natural fre- 
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quency. 
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The basic aircraft velocity response is then 
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At high elevator frequencies the velocity perturbation is essenti- 
ally negligible so that the basic aircraft velocity response is 
essentially that of a slow, second-order system of natural fre- 
quency 03^ =Cj^/y/2V damping ratio ?p= - /Cj^/27 and steady state 

gain: ^ 
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= 10.1 per rad. 


(D-41) 


s-»-0 


when full vertical gust alleviation is implemented (K_=l) but 

a 

without spring loading the flaps (K^=0) the velocity transfer 
function reduces to: 
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The zero is now shifted close to the origin of the complex fre- 
quency plane and the derivative term can no longer be neglected. 
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The essential difference however is the strong attenuation in 
the overall level of the velocity perturbation which is apparent 
from the frequency response plots and results in a steady state 
gain which is opposite in phase and reduced by a factor of 20, 
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The same technique can be applied to the phugoid approximation 
of the basic airplane's pitch angle transfer function which has 
two real zeros. One is large at approximately: 
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and one small at: 




- K 


C ^94 

Li 


0 yC 


m. 


-.06 rad/s 


At low elevator frequencies, this expression can be reduced to: 
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The steady state gain K is given by: 
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The frequency response shows a large resonant peak at the 
phugoid so that the pitch attitude response to a step elevetor 
input is dominated by the slow, lightly damped oscillations of 
large amplitude and overshoot characteristic of the phugoid mode. 
A significant change in the attitude response to low frequency 
elevator deflections occurs when full angle of attack gust alle- 
viation is introduced. The pitch angle transfer function becomes 
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And the dc gain is now: 
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To a lesser degree of accuracy this gain can be approximated by: 

C 

in ^ 

= -8.67 (D-48) 

m 

a 

The complex pair of zeros in the attitude transfer function have 
therefore approximately the same natural frequency but a smaller 
damping ratio than the characteristic phugoid roots. This imper- 
fect pole - zero cancellation is the reason for the singularity 
observed on the pitch angle Bode plot at the phugoid frequency. 
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As can be checked on the results of the analog simulation there 
is very little residual phugoid oscillation in the transient 
attitude response to a step elevator input while the steady-state 
'Value of the pitch angle is more than doubled. 

To summarize some of the most significant results of the 
preceeding analysis, the addition of the flap-auxiliary wing system, 
designed to alleviate vertical gusts, 

- Reduces the natural frequency and the damping ratio of 
the short period oscillations, 

- Considerably increases the angle of attack and pitch 
sensitivities to elevator deflections, 

- Makes elevator control of airspeed ineffective, 

- Virtually eliminates transient phugoid oscillations, 

- Produces considerable attenuation of the flight path 
angle response in the short period frequency band thus 
losing short term elevator control of altitude. 


115 




APPENDIX E 


Optimal Alleviation System Characteristics 

Complete alleviation implies no perturbation of the aircraft's 
equilibrium state. If the coefficients of the gust forcing func-: 
tions are eliminated, these conditions are satisfied. In this . 
analysis, the required relations are developed without regard to 
the dyneimic behavior of the resulting configuration. To obtain 
these relations, vertical and horizontal gust alleviation can be, 

considered independently. Since the former is more important, it 

■ > . ^ ' 

is treated first. 

If an aircraft makes use of control surfaces to eliminate 
the gust forcing terms in the three equations of motion, these 
terms and those introduced by the control surfaces must offset one 
another. In particular, if wing flaps are employed, the following 
relations must hold; 

A-6 + A a„ = 0 (E-1) 

6 Ug g 

where the elements of the three-dimensional A vectors are given by 
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Except for drag, each of the forcing terms involves the gust-induCed 

angle of attack and its derivative. 

This set of equations is overdescribed to yield a solution for 

\ 

flap deflection as a function of gusts alone, and yields a set of 

conditions on the elements of the A vectors which also must be met 

to satisfy equation (E-1) if a is to remain arbitrary. For example, 

if the drag equation is solved for 6 as a function of a and this 
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solution is inserted in the lift equation, the relationship of the 
A "vector elements can be expressed by the determinant given in 
equation (E-2)i: Similarly, if the lift and pitching moment equations 
are treated in this manner, equation (E-3) results. 
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The second of these determines the absorption of vertical gusts in 
the short period frequency band and as such has a major impact on 
the overall smoothness of flight in rough air. 

If gust alleviation is to be achieved independent of frequency, 
it is seen that the required conditions can be satisfied by setting 
each of the resulting polynomial coefficients to zero. This is 
equivalent to satisfying the following three relations. 
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Condition (E- 4) is automatically satisfied if the contributions of 
the wing-fuselage combination to the rate derivatives are negli- 
gil?le relative to those of the tail. With this assumption the 
stability derivatives take the following forms; 
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where is the lift curve slope of the tail, nondimensionalized 

Ij 

“t 

with respect to wing area and dynamic pressure. The moment deri- 
vatives have the same structure. 

Conditions (E-5) and (E-6) reduce to 
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These relations were originally developed by Phillips and Kraft^. 

The lift portion of Equation (E-l)is solved for the flap deflec- 
tion transfer function: 
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The definition of the stability derivatives and condition (E-8) can 
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be used to show that: 




= 0 
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with the consequence that flap deflection becomes independent of 
frequency. 



(E-10) 


When the conditions for optimality are satisfied, the lift and pitch- 
ing moment of both wing and tail become independent of angle of 
attack. The flaps will then act as perfect acceleration alleviators 
in the short period frequency range, in the sense that operation of 
the flaps to offset the accelerations produced by high frequency 
vertical gusts will result in no pitching motions. Furthermore, 
the flap deflection required under these conditions can be imple- 
mented by a simple passive aeromechanical system of high natural 
frequency so that the flaps respond fast compared to gust perturba- 
tions . 

Condition (E-2) determines the alleviation system characteristics 
necessary to filter out completely the low frequency component of 
vertical gusts. Elimination of the polynomial coefficient yields 
the following two relations: 
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These are equivalent for the system considered because of 
Equation (E-9) and both reduce to relation (E-13) which defines the 
optimal drag characteristics of the lift control surfaces: 



(E-13) 


The conditions (E-7) , (E-8) and (E-13) for optimal vertical gust 
alleviation can be rewritten in simple form with the aid of the 
total derivative. 
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Examination of these relations shows that the chordwise distri 
bution of incremental lift associated with control surface deflec- 
tion should produce zero pitching moment about the aerodynamic cen-\: 
ter of the wing. This is impossible to realize in practice although 
coupling of the elevator with the flaps could be used to offset this 
natural pitching tendency. Also, in order for the angle of attack 
at the tail to be independent of that of the wing as required by 
the optimal downwash condition, the flaps should produce downwash 
on the stabilizer when sensing an upward gust which indicates that 
the spanwise distribution of the incremental lift should be concen- 
trated towards the wing tips. Effectively this means that the lift 
control surfaces should be of the aileron type i.e. located outboard 
The drag condition implies the need for employing a device to 

4 

counteract . Phillips has recognized the impact of this para- 

meter for a S^OL aircraft. However, this aircraft exhibits much 
larger values for the lift curve slope and incremental drag due to 
flap deflection than does the Cessna 172. As a result the ratio 
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of drag and lift increments due to a vertical gust is more than . 
an order of magnitude smaller for the Cessna 172 than for the 
STOL aircraft and longitudinal accelerations introduced by the 
gust alleviation systems should be tolerated easily in the light 
aircraft. 

A similar analysis can be performed with respect to hori- 
zontal gusts, yielding another set of optimal flap characteristics. 
Stability derivatives with respect to rate of change of gust velo- 
city are introduced to measure the effects of the lag in dynamic 
pressure on the tail when the aircraft penetrates a region of 
horizontal gust gradients: 

Ct =-2i Cr (E-15) 

. . "g 

where is the lift coefficient of the tail normalized with re- 
"t . 

spect to wing area and dynamic pressure. is the corresponding 

moment coefficient derivative about the center of gravity. Ideally 
the gust alleviation system should completely absorb the incremental 
forces and moments produced by the horizontal component of atmos- 
pheric turbulence in order to reduce the net forcing terms in the 
equation of motion to zero over the relevant gust frequency 
spectrum. 
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where 
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The imposed conditions on the coefficients of the A vectors can be 
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expressed by. 
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The analysis developed previously in the case of vertical gusts 
can be followed step by step using the S 2 ume assumptions concern- 
ing the rate derivative and yields the following three optimality 
conditions , 
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where the total derivatives are defined by: 
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In this case, the optimal flap deflection also is found to be inde- 
pendent of gust frequency. 
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This response can be achieved physically by loading the flaps with 
a spring in such a way that the added hinge moment does not vary 
with flap deflection or dynamic pressure. 

The first of the optimality conditions might be satisfied by 
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a4 jutting the loading on the tail so that: 

(E-21) 
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w^ich generally is positive. The third condition specifies the 
drag characteristics of the flaps required to eliminate the longi- 
tudinal accelerations produced by low frequency horizontal gust, 
and could in principle be satisfied through the use of spoilers 
coupled to the flaps. 

A close examination of the conditions for optimal horizontal 
and vertical gust alleviation shows them to be incompatible. The 
zero flap pitching moment requirement for optimal horizontal gust 
alleviation cannot be reconciled with the zero effective pitch 
stiffness condition for optimal vertical gust alleviation unless 
the basic aircraft is neutrally stable. Furthermore optimal verti- 
cal gust alleviation requires that flaps produce upwash on the tail 
when deflected downward whereas they should produce downwash for 
optimal horizontal gust alleviation. 

Since vertical gust alleviation is the more useful, the 
stability characteristics of the corresponding optimal configuration 
are examined. It is shown that if the flap-auxiliary wing system 
is designed for optimal vertical gust alleviation, as defined in 
the previous section, the augmented aircraft is marginally stable, 
and the short period mode of oscillation degenerates into a pair 
of real roots, one of which is at the origin of the complex fre^- 
quency plane. Marginal stability is implied directly by one of. the 
Ojptimality conditions previously developed, 
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When all of the optimality conditions for vertical gust 
alleviation are met, and horizontal gust alleviation is neglected, 
the characteristic determinant reduces to: 
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Expansion of the determinant neglecting small terms yields the 
fourth degree polynomial 


P(s) = 4y^I s^ - 8]X^C^ + 4yCjjj s^ - 2C s (E-25) 

^ q q u q 

The rigid body modes of the augmented aircraft are clearly separ- 
ated so that this characteristic equation can be factored approxi- 
mately into: 
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The short period mode is thus degenerated into a pair of real roots, 
one of which is at the origin of the complex frequency plane (s=0) 

m 


while the other is large and stable at s = 


-6.30 rad/s. 


As expected, the phugoid is relatively unaffected by the in- 
troduction of the angle of attack compensation system. The na- 
tural frequency of the mode is reduced by 30% =Cj^/2y:.166 rad/s) 

P 

while the damping ratio is increased, (C ="C /2C_=0.14) so that 

P *u u 

the time to halve amplitude remains unchanged. 
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The degenerate short period mode associated with the loss 
of pitch stiffness resulting from the requirements for optimal 
vertical gust alleviation is obviously unacceptable from the 
point of view of stability. It should be emphasized here that 
applying gust alleviation conditions is contrary to normal stability 
requirements. In any practical system a satisfactory compromize 
must be reached. 
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